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Stimulated Raman scattering (SRS) in ultra-high-Q surface-tension-induced spherical and chip- 
based toroid microcavities is considered both theoretically and experimentally. These microcavities 
are fabricated from silica, exhibit small mode volume (typically 1000 /^m^) and possess whispering- 
gallery type modes with long photon storage times (in the range of 100 ns), significantly reducing 
the threshold for stimulated nonlinear optical phenomena. Oscillation threshold levels of less than 
100 /I- Watts of launched fiber pump power, in microcavities with quality factors of 100 million are 
observed. Using a steady state analysis of the coupled-mode equations for the pump and Raman 
whispering-gallery modes, the threshold, efficiencies and cascading properties of SRS in UHQ devices 
are derived. The results are experimentally confirmed in the telecommunication band (155Qnm) 
using tapered optical fibers as highly efficient waveguide coupling elements for both pumping and 
signal extraction. The device performance dependence on coupling, quality factor and modal volume 
are measured and found to be in good agreement with theory. This includes analysis of the threshold 
and efficiency for cascaded Raman scattering. The side-by-side study of nonlinear oscillation in 
both spherical microcavities and toroid microcavities on-a-chip also allows for comparison of their 
properties. In addition to the benefits of a wafer-scale geometry, including integration with optical, 
electrical or mechanical functionality, microtoroids on-a-chip exhibit single mode Raman oscillation 
over a wide range of pump powers. 



I. INTRODUCTION 



Ultra-high-Q ( UHQ) surface-tension-induced micro- 
cavities (which we subsequently refer to as STIMs) com- 
bine small modal volume with some of the highest op- 
tical quality- factors (Q) to date of nearly 10 billion [1], 
and are of interest for a variety of studies ranging from 
fundamental physics such as cavity quantum electrody- 
namics [2] [3] [4] to applied areas such as low threshold 
and narrow linewidth lasers [5] [6] [7] , nonlinear optical 
oscillators [8] [9] , as well as high-sensitivity transducers for 
biochemical sensing[10]. These silica microcavities fea- 
ture whispering gallery type modes (WGMs) and rely 
upon exquisite smoothness at the cavity dielectric bound- 
ary to attain ultra-high-Q performance that is typically 
in excess of 100 million. For nonlinear optical studies, 
strong resonant buildup of energy in niicroscale volumes 
significantly reduces the threshold for nonlinear optical 
effects to occur. This was recognized in the pioneering 
work of Chang [11] [12] and CampiUo [13] [14] [15] who ob- 
served and studied a variety of nonlinear optical effects in 
ultra-high-Q liquid microdoplets. Their work used free- 
space illumination to optically pump the microdroplets 
and thereby induce Raman oscillation [12] [13] [15], cas- 
caded Raman scattering [12] and Brillouin scattering[16]. 
Silica UHQ STIMs provide a far more stable and robust 
microcavity in comparison with liquid microdroplets. 
However, despite numerous studies on these devices over 
the past decade [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] the 
observation of nonlinear phenomena (beyond thermal ef- 
fects) in these devices, had been limited to one report 
on Kerr-induced wavelength shifts at low temperatures 



[20]. The advent of pumping and signal collection using 
fiber taper couphng methods [27] [28] [29] [30] proved an 
important turning point in access to nonlinear phenom- 
ena in this important micro-cavity system. Fiber tapers 
provide remarkably efficient coupling to and from the 
UHQ silica sphere system [28] [30] and enable direct ac- 
cess to the technologically important optical fiber trans- 
port medium. The measured fiber-coupled threshold for 
a variety of nonlinear phenomena in taper-coupled, silica 
microspheres are lower than for any other nonlinear oscil- 
lator reported to date. Silica microsphere Raman lasers 
with ultra-low threshold levels of only 62 /i- Watts [31] 
have been demonstrated. Compared to microdroplets 
these devices allow stable and long term observation of 
nonlinear optical effects in microcavities. Cascaded Ra- 
man lasing in these devices of up to 5 orders has also 
been observed [32] . The tapered optical fiber in these 
experiments functions to both pump WGMs as well as 
to extract the nonlinear Raman fields. In addition, the 
tapered-fiber coupling junction is highly ideal[28], mak- 
ing it possible to strongly overcouple ultra-high-Q cav- 
ities with negligible junction loss. This feature allows 
for the observation of very high internal differential pho- 
ton conversion efficiencies approaching unity. Whereas 
microspheres are both compact and efficient nonlinear 
oscillators, their fabrication properties lack the control 
and parallelism typical of microfabrication techniques. 
Recently-developed ultra-high-Q toroid microcavities on- 
a-chip [33] provide a UHQ silica device with performance 
equivalent to a microsphere. UHQ toroids have several 
advantages over spheres including being wafer-scale de- 
vices that can be fabricated in parallel as dense arrays 
or integrated with electronics or other optical functional- 
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ity. In this paper, we demonstrate and analyze nonlinear 
Raman oscillation in both microsphere and microtoroid 
on-a-chip structures. In addition to studying Raman os- 
cillation in microcavities in detail, this work allows us 
to compare the performance and properties of toroidal 
and spherical microcavities. In addition to their fabri- 
cation and integration advantages, it will be seen that 
microtoroids also have performance advantages in com- 
parison to microspheres. This includes a reduced number 
of supported azimuthal modes, which allows observation 
of single-mode Raman oscillation over a large range of 
pump powers [34] (of critical importance in practical ap- 
plications), as well lower Raman threshold due to a re- 
duced mode volume compared to a spherical cavity. 

The paper is organized as follows. The first section 
will give a brief introduction into the fabrication, cou- 
pling and optical properties of ultra-high-Q toroid on-a- 
chip and spherical microcavities. The second section will 
present a model for Raman lasing in a waveguide-coupled 
whispering-gallery-mode resonator. There, we derive the 
expression for the threshold and the efhciency of the con- 
version process. In the ultra-high-Q regime resonances 
are easily split into doublets due to intermode coupling 
of the degenerate clockwise and counterclockwise prop- 
agating whispering gallery modes [35]. The effect of this 
intermode coupling on stimulated Raman scattering is 
also considered in this section. The case of cascaded Ra- 
man oscillation in which Raman signals serve to pump 
and generate higher-order Raman waves is also treated in 
the section. In the fourth section, the theoretical results 
are compared with experimental studies of the depen- 
dence of Raman threshold on quality- factor, mode vol- 
ume and waveguide loading. Experimental results con- 
cerning power and efficiency during cascaded operation 
are presented and compared to theory and found to be 
in good agreement with theoretical predictions. Finally, 
a comparison of microtoroid oscillation properties with 
those of microspheres is presented in the fifth section. 



II. ULTRA-HIGH-Q 
SURFACE-TENSION-INDUCED 
MICROCAVITIES 



Surface-tension-induced microcavities such as micro- 
droplets, exhibit a superb, cavity surface finish (typically 
nanometer surface roughness [1]) leading to whispering 
gallery type modes with some of the highest optical qual- 
ity factors recorded to date. In the work presented here, 
both spherical and toroid microcavities on-a-chip made 
from silica are investigated. Both types of structures ex- 
hibit Q-factors in excess of 100 million. Briefly, the fabri- 
cation of both micro-sphere and micro-toroids relies upon 
surface tension to induce collapse of a given silica preform 
into the final cavity shape. Because the cavity fabrica- 
tion involves a temporary liquid state the surface finish of 




FIG. 1: Optical micrograph of a spherical (left) and toroid 
microcavity (right) coupled to a tapered optical fiber. 



the final "solid" microcavity is excellent. In the case of a 
microsphere, the preform used here is an optical fiber tip 
which is heated and melted with a carbon-dioxide laser 
(10.6 /im wavelength). Surface tension causes the silica 
fiber tip to contract and form a spherical microcavity, 
while the remainder of the fiber stem serves as a holder 
for sphere-positioning. Light within the sphere is con- 
fined near an equatorial plane by continuous total inter- 
nal reflection at the cavity interface. Microtoroids, on the 
other hand, use wafer microfabrication techniques involv- 
ing a combination of lithography and etching combined 
with a final selective reflow process using a CO2 laser. 
The details of this fabrication process are reported in 
reference [33] . In the case of toroid microcavities the pre- 
form consists of a microfabricated silica disk supported 
by a silicon pillar. Illumination of the disk using a CO2 
laser induces selective reflow of the silica. Surface ten- 
sion causes the disk preform to collapse into a toroidal 
periphery, thereby creating the resonant cavity. Figure 
1 contains an optical micrograph of a microsphere and a 
microtoroid on-a-chip. 

Highly efflcient evanescent coupling to these micro- 
cavities can be achieved by use of tapered optical fibers 
[30] [28], which are fabricated by melting and adiabati- 
cally tapering a standard, telecommunication fiber un- 
til a waist diameter of approximately 1 to 2 microns is 
reached. When brought into proximity of the WGM res- 
onator and when the taper waist diameter is chosen to 
phase match to the WGMs, efficient coupling both to 
and from the microcavity can be achieved. Coupling 
can be described by the normalized coupling parameter 

Kp — (7^) which describes the ratio of intrinsic res- 
onator lifetime tq to the external (coupling related) life- 
time Tex- Following the standard conventions [36], under- 
coupling is denoted by K < 1, overcoupling by K > 1, 
and critical coupling (the point of vanishing waveguide 
transmission) is denoted by K = 1. The ability to achieve 
strong overcoupling is important as it exemplifies the 
very ideal nature [28] of the taper-microcavity, coupling 
junction. In addition to efficient excitation and extrac- 
tion of optical power, the fiber taper provides direct cou- 
pling to optical fiber, thereby further facilitating labora- 
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FIG. 2: Cavity rindown measurement of a 50- micron-diameter 
toroid microcavity coupled to a tapered optical fiber. At t=0 
the laser is gated off, and (after an 8 ns fall-time) the trans- 
mission is entirely due to the cavity decay field. Inset: A 
logarithmic plot to infer the decay time yields 75 ns at the 
critical point. 



of 100 MHZ, makes overlap of cavity modes with the 
Brillouin gain spectrum unlikely and was not observed 
in this work. In the case of parametric oscillation me- 
diated by the Kerr-nonlinearity (i.e. four- wave mixing), 
energy conservation requires a triple resonance condition 
for signal, idler and pump mode, which can be satis- 
fied only under certain conditions [37]. These additional 
constraints for oscillation based on the Kerr nonlinearity 
or Brillouin gain, cause stimulated Raman scattering - 
which is intrinsically phase-matched- to be the dominant 
microcavity nonlinear optical effect. For simplification, 
we assume that the pump wavelength and the Raman 
wave are on resonance and use the slowly varying enve- 
lope approximation. 

(1) 



tory measurements. In the case of microtoroids on-a-chip 
the taper is also crucial as a means of " probing" the whis- 
pering gallery devices which are within a few microns of 
the silicon wafer surface. 

Using tapered optical fibers, resonator quality factor 
can be inferred from either linewidth measurements or 
cavity ringdown experiments [33]. Figure 2 shows a cav- 
ity ringdown measurement on a 45-/^r7i-diameter micro- 
toroid at 1550 nm, exhibiting a critically coupled Q of 
10®. In this measurement the WGM was excited on reso- 
nance, and the fiber-taper adjusted to the critical point. 
The cavity lifetime can then be inferred by gating-off the 
excitation laser, and recording the cavity decay signal. 
The critically coupled Q (including the waveguide cou- 
pling contributions) was 100 million, and when correcting 
the Q for the waveguide loading, as described in section 
III.B, an intrinsic quality factor of 3.7 x 10® is inferred. 



Here E signifies the slowly-varying amplitude of the 
pump and Raman WGM modes of the cavity and s de- 
notes the input wave. The excitation frequency of the 
pump mode and resonant Raman mode is given by ojr 
and Wpand r is the total lifetime of photons in the res- 
onator, which is related to the quality factor hy Q — uj ■ t. 
The coupling coefficient k denotes the coupling of the 
input pump wave s to the cavity whispering-gallery- 
mode The relation k= i/— associates the cou- 

pling coefficient with a corresponding lifetime, such that 
- — — — I- — [36] . Since the Raman effect will excite both 

T Tex TO L J 

eigenmodes of the cavity (co- and counter-directionally 
propagating modes), equal amplitude emission occurs 

along both waveguide directions given by |sr|^ = 2t • 
The Raman intra-cavity gain coefficient is denoted as g^j, 
which can be related to the more commonly used gain co- 
efficient (7fl,(measured in units of m/Watt) by. 



III. THEORETICAL ANALYSIS OF RAMAN 
SCATTERING IN HIGH-Q MICROCAVITES 

A. 1st Order Raman Scattering in microcavities 

Raman scattering in a waveguide coupled microcav- 
ity can be described classically by using coupled mode 
equations for the pump and Raman fields with non- 
linear Raman coupling terms. Other nonlinear effects, 
which can compete with Raman scattering, such as Four- 
wave-mixing or Brillouin scattering are not considered 
in this analysis as the microcavity poses stringent fre- 
quency matching constraints on these processes making 
their observation difficult. In the case of Stimulated Bril- 
louin Scattering, the narrow gain bandwidth in the range 
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where Veff is the effective modal volume [38], and E is 
the electric field vector [2] . The effective mode volume ac- 
counts for the intensity dependent gain, and for silica mi- 
crospheres and microtoroids has approximately twice the 
value than the energy related definition of mode volume. 
Steady state analysis of the coupled mode equations, re- 
sults in a clamped, cavity, pump field above threshold. 
This clamping alters the coupling of pump power to the 
resonator, and, in turn, the pump power dependence of 
Raman laser power such that the following square root 
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dependence results. 
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The physical origin of this square root dependence of the 

pump- to- Raman conversion can be viewed as a "pump- 
ing inefficiency" i.e. the coupled pump power does not 
increase linearly with launched fiber power. The nonlin- 
ear dependence of coupled pump power can be illustrated 
by considering a pump wave that is initially critically 
coupled to the resonator. As noted above, critical cou- 
pling features complete transfer and dissipation of power 
from the resonator (i.e., zero transmission). In terms 
of the fields involved in coupling both to and from the 
resonator, critical coupling results from the destructive 
interference of the cavity leakage field with the trans- 
mitted, pump field (i.e., the portion that does not cou- 
ple to the resonator from the waveguide). Once the on- 
set of Raman lasing is reached, the cavity pump field is 
clamped at the threshold value resulting in a fixed cav- 
ity pump leakage field. Subsequent increase in launched 
pump power will imbalance the leakage and the transmit- 
ted pump fields, giving rise to finite transmission and a 
shift away from the critical point. The pump coupling to 
the resonator is thereby less and less efficient as the pump 
field is increased. The expression for the Raman thresh- 
old pump power can be factorized into terms involving 
modal volume, waveguide-cavity coupling strength and 
cavity lifetime (or Quality factor). To facilitate separa- 
tion of the coupling and intrinsic lifetime dependence, 
we use the dimensionless normalized coupling parameter 

Kp = (^7^) • III the ideal case of a single- mode waveg- 
uide coupled to a whispering-gallery-mode the waveg- 
uide transmission as a function of coupling is given by 

T = (tT^) ^^'^ ^ typically varies exponentially with 
the waveguide- microcavity "coupling gap" distance [36]. 
Using these definitions and under the assumption of equal 
coupling properties and photon lifetimes for both the 
pump and Raman mode, i.e. Kp = Kr = K and tr = Tp 
,the threshold expression is given by: 
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(4) 



Here we have also introduced C{T) which is a possible 
correction factor to account for intermode coupling of the 
degenerate clockwise and counterclockwise propagating 
whispering-gallery modes. This factor will be explained 
in the next section. The threshold expression follows an 
inverse square dependence on the quality factor. This re- 
fiects the fact that an increase in Q will cause a twofold 
benefit in terms of both reducing cavity round trip losses 
that must be overcome for threshold as well as increasing 
the Raman gain, due to the intensity dependence of the 
Raman gain coefficient on the pump field. In addition, 
the equation shows that the threshold scales linearly with 
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FIG. 3: External and internal pump to Raman differential 
conversion efficiency. The external efficiency refers to conver- 
sion of launched pump power to Raman output power, while 
the internal efficiency is the coupled pump power to Raman 
output power. As evident, the external efficiency can exceed 
unity in the overcoupled regime 



the modal volume. Both the coupling and mode volume 
dependence of the Raman threshold are examined exper- 
imentally in the next section. When analyzing the cou- 
pling dependence under the assumption of equal Raman 
and pump cjuality factors and coupling factors, the min- 
imum threshold occurs when Q™" = 2 • Qo or K"^™— 
i.e. in the undercoupled regime with finite waveguide 



transmission of T™'"= 
threshold is given by. 
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It is worth noting that at this coupling condition, the cir- 
culating pump-power in the resonator is not maximum. 
This can be understood since minimum threshold repre- 
sents an optimal balance of both pump coupling and Ra- 
man mode coupling loss. The conversion of pump power 
to Raman power above threshold can be characterized by 
the differential slope efficiency. The bidirectional exter- 
nal differential slope efficiency 7]^^ is derived by lineariz- 
ing the expression for Pr near the threshold condition 
and is given by. 
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Figure 3 shows the differential slope efficiency as a func- 
tion of coupling strength. It is noteworthy that it ap- 
proaches the value of 2^ in the limit of strong overcou- 
pling {-^ = oo). Surprisingly, this value exceeds unity, 

indicating that every waveguide pump photon added 
above threshold, is converted to more than one Raman 
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photon. This rcsuh can be understood by again eonsid- 
ering the nonhnear dependence of coupled pump power, 
except this time in the over-coupled regime. In par- 
ticular, the differential increase in coupled pump power 
grows more quickly in the overcoupled (more slowly in 
the undercoupled regime) than the differential increase 
in launched pump power. This leads to the interesting 
effect that the differential photon conversion efficiency 
can exceed unity. Taking into account the nonlinear de- 
pendence of coupled pump power by defining the internal 
differential efficiency jjj^t as the coupled (as distinct from 
the launched) pump-to-Raman power, the efficiency ap- 
proaches, as expected , the value ^in the limit of strong 
over coupling. 

dPRaman ^ij / 1\ ^ 

Vint = -T5 = 7^\^ + l?) 

U-i coupled '^p \ ^ / 

Figure 3 shows both the internal and external differential 
Raman conversion efficiencies as a function of coupling 
strength. 
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FIG. 4: Power correction factor (which is the ratio defined 

by maximum circulating power in the absence and presence 
of modal coupling) as a function of F-parameter. The inset 
shows a frequency scan of a doublet of an ultrar-high-Q mode. 



B. The effect of intermode coupling on stimulated 
Raman scattering 

In the ultra- high-Q regime, resonances of a whispering- 
gallery-type microcavity are often split into doublets[17] 
[35]. This splitting is due to coupling of the degenerate 
clockwise and counterclockwise propagating modes by ei- 
ther intrinsic or surface scattering centers. The modified 
coupling properties have been extensively studied [39]. 
Here we briefly summarize the results of this study and 
analyze the effect of intermode coupling on Raman scat- 
tering. The extent to which intermode coupling modifies 
the waveguide coupling properties can be described by 
the dimensionless intermode coupling parameter F = ^ , 
where 1 /j is the rate of coupling of the degenerate clock- 
wise and counterclockwise modes. It is easily measured 
as the linewidth normalized splitting in the undercoupled 
regime (see figure 4 inset). The presence of intermode 
coupling has several consequences. First, the critical 
point (as defined by vanishing waveguide transmission) 
is shifted towards a coupling point, which, using the con- 
ventional coupling terminology, is considered overcoupled 
as Tex < To. The point of vanishing transmission occurs 
at: 



Third, the leakage of the cavity field into the back- 
wards waveguide direction causes a reduction of the cav- 
ity buildup factor with respect to the ideal case in the ab- 
sence of intermode coupling. In the limit of strong modal 
coupling, the cavity buildup factor is reduced by a factor 
of 2, which subsequently causes a twofold increase in the 
threshold necessary to achieve Raman lasing. Figure 4 
shows the circulating power correction factor C{T) as a 
function of the dimensionless intermode coupling param- 
eter F. In the presence of modal coupling the waveguide 
coupling condition for minimum threshold experiences a 
slight shift towards overcoupling with the maximum shift 
occurring at F « 1.52. In the regime of very strong modal 
coupling the condition of minimum Raman threshold ap- 
proaches again the original condition K = 1/2. 



Kcrit = + 1 



(8) 



Second, the shifted critical point is accompanied by a 
maximum refiection into the backwards direction of the 
waveguide (i.e. the contra-directional waveguide mode). 
The magnitude of the refiection at the above modified 
critical point is given by 



(9) 



Analysis of Cascaded Raman Scattering in 
high-Q microcavities 



The first Raman field can itself act as a secondary 
pump field and generate further Raman modes. This 
process of cascaded Raman scattering can be described 
by including higher order coupling terms into the cou- 
pled mode equations of pump and Raman fields as shown 
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below [32]. 
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where A'^ is the Raman order. To find the correspond- 
ing thresholds and output powers for these higher order 
processes, the set of equations can be solved iteratively in 
steady state. Here, we have introduced the dimensionless 
coefficients &. 
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The general solutions for the threshold of the even and 
odd order TV*'' Raman modes are given by the following 
expressions. As in the previous section, we have assumed 
equal coupling strengths and intrinsic Q factors for the 
pump and Raman modes. 
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As evident from these expressions even and odd order 
stokes fields exhibit different threshold powers as a func- 
tion of stokes order (N). When considering Raman scat- 
tering in silica at optical frequencies, one can approxi- 
mate the above expressions by taking Ci ~ 1 since the 
Raman shift is small compared to the optical frequency. 
In addition it is assumed that the mode volume is wave- 
length independent. Under this assumption, the thresh- 
old expressions reduce to: 
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It follows that the threshold for cascaded Raman oscil- 
lation exhibits a cubic dependence on Raman order N. 
The emission power dependences vary depending upon 
whether the highest order wave is even or odd. For the 
odd order case, all odd orders increase as the square root 
of the pump power and even orders are clamped. For 
the even order case, all even order lines increase linearly 
with pump power while odd orders are clamped. Figure 
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FIG. 5: Theoretical plot of cascaded Raman output power 
versus pump power for up to fourth order Stokes power. Even 
and odd ordered Raman modes exhibit alternating clamping 
behaviour. 



5 illustrates this behavior showing the Raman output for 
several stokes orders as a function of input pump power. 
The analytic expressions for the Raman output power in 
these cases are given by: 
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The differential power conversion efficiencies can be ob- 
tained by linearizing the above expressions near the 
threshold condition. The external and internal differ- 
ential efficiencies decrease steadily as a function of stokes 
order {N). For optical frequencies that are much larger 
than the Raman shift, the external differential conversion 
efficiency reduces to: 
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For high order (N) Raman fields, the external differential 
conversion efficiency thus follows a 1/iV^-dependence. 



IV. EXPERIMENTAL STUDY OF 
STIMULATED RAMAN SCATTERING IN 
ULTRA-HIGH-Q MICROSPHERES 

We have observed stimulated Raman scattering us- 
ing fiber-taper-coupled, ultra-high-Q, surface-tension- 
induced silica microspheres. Tapered optical fibers pro- 
vide both a very efficient and practical means for pump- 
ing and laser signal extraction through the same fiber. 
The experimental setup in this study is identical to the 
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FIG. 6: Single mode Raman lasing in a UHQ spherical mi- 
crocavity. The pump wavelength is located at 1550 nm and 
Raman lasing appears at 1660 nm. The inset shows the Ra- 
man frequency shift (designated by the dotted line) overlaid 
with the Silica Raman gain spectrum. 
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FIG. 7: Spectrally resolved Raman emission of the ultra-high- 
Q microsphere shown in the previous figure, revealing that 5 
eccentricity split azimuthal modes were oscillating simultane- 
ously. The pump power in this measurement was adjusted 
to slightly above Raman oscillation threshold. The assym- 
metric line-shape of the Raman modes is caused by repeat- 
edly scanning the pump laser through the pump resonance, 
which causes a redshift of the cavity resonances, due to heat- 
ing and subsequent expansion of the cavity (thermal bistabil- 
ity as noted in [17]). 



one reported earlier [31] and the tapered fiber was at- 
tached to a piezoelectric stage which allowed precise con- 
trol of the taper-microcavity coupUng gap (20 nm resolu- 
tion) and variation of the waveguide-resonator coupling 
strength. The Q-factors obtained in the microspheres 
were typically in the range of 1 — 2 x 10* and taper in- 
sertion loss (fiber to fiber) was usually less than 5%. 

To excite the UHQ modes, we used a narrow-linewidth, 



external-cavity laser emitting in the 1550 nm band. In 
the case of microspheres, the WGM field spatial forms are 
known analytically and characterized by the radial, angu- 
lar, azimuthal and polarization mode number (n, I, m,p). 
Due to fabrication- induced eccentricity, the (2^-1- l)-fold- 
degeneracy of the azimuthal modes is lifted yielding a 
complex mode spectrum. For the spheres considered in 
this work the eccentricity-induced splitting was in the 
range of several GigaHcrtz. Microtoroids on-a-chip, on 
the other hand, possess a significantly reduced mode 
spectrum due to the strong azimuthal modal confine- 
ment provided by the toroid geometry. This both simpli- 
fies their spectra and enables operation of microtoroid- 
Raman lasers in the desirable single mode regime, which 
is of significant practical importance. 

Stimulated Raman oscillation was observed by pump- 
ing a single WGM and monitoring the transmission us- 
ing an optical spectrum analyzer. Once the threshold 
for SRS was exceeded, lasing modes in the 1650-nm 
band could be observed, in correspondence with the peak 
Raman gain which occurs downshifted in frequency by 
approximately 14 THz relative to the pump frequency 
(wavelength shift of approximately 110 nm). Figure 6 
shows Raman emission for an ultra-high-Q microsphere. 
The Raman emission with respect to the gain peak is pro- 
vided in the inset of figure 6. Since the fundamental whis- 
pering gallery modes (n = 1,£ = m,p = TM) are most 
tightly confined (i.e., smallest mode volume), Raman las- 
ing is expected to occur first for these modes. The pres- 
ence of nearly degenerate azimuthal modes in a spherical 
microcavity (i.e., weak eccentricity splitting), causes si- 
multaneous oscillation on several azimuthal modes. Fig- 
ure 7 shows a higher resolution spectral scan of the spec- 
trum in figure 6. Several azimuthal modes can be ob- 
served to be oscillating simultaneously. 

The threshold formula predicts a strong dependence 
of the Raman threshold on waveguide coupling. Figure 
8 shows the measured Raman threshold as a function of 
taper-microcavity gap distance for a fundamental WGM. 
The solid line is a fit using the theoretical results from the 
first section (equation 3) and assuming equal pump and 
Raman Q factors and coupling conditions. The minimum 
threshold does indeed occur undcrcoupled with finite ta- 
per transmission at T = 13% in good agreement with the 
theoretically predicted value of T = 1 /9. 

As a further verification of the threshold formula we 
compared the theoretical minimum threshold value with 
the observed value. The quality factor and the mode 
splitting of the whispering-gallery mode were measured 
by performing a linewidth sweep in the undercoupled 
regime, where the backscattering-induced doublet struc- 
ture is most pronounced. These measurements yielded 
Qo = Ix 10* and F = 2 . The size of the microsphere was 
inferred from the free spectral range, AA = 10.5 nm i.e. 
50 /xm diameter, where the free-spectral range denotes 
here modes with successive angular mode number £. The 
mode volume was calculated using analytic expressions 
based on estimated mode numbers for the fundamental 
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FIG. 8: The dependence of Raman threshold on loading as 
evidenced by the taper-micro-cavity gap distance. The taper- 
microcavity gap is measured in microns with respect to the 
critical point (defined as vanishing taper transmission), and 
positive separation corresponds to overcoupling. The mini- 
mum threshold is 62 /iW of launched pump power (measured 
before the taper which includes the insertion loss of the op- 
tical fiber taper). Minimum threshold occurs ca. 120 nm 
undercoupled (corresponding to a transmission of 13%). 




microsphere Size (|im) 

FIG. 9: The Raman oscillation pump threshold of a spherical 
UHQ microcavity as a function microsphere diameter. To 
compare microspheres with different resonance characteris- 
tics, the Raman threshold was Q-normalized to Q = 1.0 x 10* 
and modal coupling normalized to F = 0. The inset shows 
a double logarithmic plot of the data. A linear fit obtains 
a dependence of R^'^^ which is in close agreement with the 
theoretically expected dependence of R^'^^. 



WGM (n = l,p = TM, m^i,l^ 139). Calculations [40] 
[41] yielded a modal volume of ca.1300 //m^. The over- 
lap factor for pump and Raman mode is assumed to be 
unity. Using these values, the theoretically expected min- 
imum threshold is given by 50 fiWatts which is in good 
agreement with the experimentally measured value of 62 
fiWatts. 

The dependence of the Raman threshold on the modal 
volume was also investigated. For this measurement the 
Raman threshold was measured for micro-sphere res- 
onators having varying diameters in the range of ca. 
25-120 /iTO. (For diameters smaller than 25-/im ther- 
mally induced wavelength shifts lead to pumping insta- 
bilities.) In this diameter range, the mode volume follows 
an approximately quadratic dependence (actual inferred 
exponent is cx R^'^^ [41] [40]) on the sphere radius. 
As an aside, for smaller spheres the mode volume de- 
viates from this behavior and ultimately, for very small 
diameters, the mode volume increases due to weaken- 
ing of the whispering gallery confinement [40] . The min- 
imum mode volume occurs for a radius of 6.9/im (for 
£ = m = 34) [40] for 1550 nm wavelength and the mode 
volume is Vmin = 173.1/xm'^. However, this size is not op- 
timum for stimulated Raman scattering as the additional 
benefit of reduced mode volume is more than offset by 
the significant decrease in Q factor to 10^ [41] (Threshold 
power (X V/Q'^). For small mode volumes, it has been 
predicted that the gain coefficient can exhibit a depen- 
dence on mode volume due to cavity QED effects [14]. 
However, in the case of stimulated Raman Scattering in 
silica microspheres these effects are not expected to be 



observable [42] and an approximately quadratic relation- 
ship due to the mode volume is predicted. 

Figure 9 shows the experimental results for thresh- 
old versus microsphere diameter. In this experiment, 
the minimum Raman threshold P™'", the microsphere 
size (as inferred from the free-spectral-range), the intrin- 
sic Q in the pump band (Qq) and the intermode cou- 
pling parameter (F) were measured. To extract the volu- 
metric dependence of the Raman threshold using data 
from cavities having different resonant characteristics, 
the threshold data were normalized to the set of param- 
eters (Qo = 10^^ = 0,gR = 5^^^^^). The result of this 
procedure is shown in figure 9. The data indeed show 
a quadratic dependence on R (the actual fitted exponen- 
tial from a double logarithmic plot is 1.95 and is in good 
agreement with the expected value of 1.83) and confirm 
the linear relationship of the Raman threshold on the 
mode volume as predicted by equation 4. 

A. Cascaded Raman Scattering in ultra-high-Q 
microspheres 



In addition to first order Raman scattering, we have 
also observed cascaded Raman scattering. Figure 10 
shows a typical cascaded Raman spectrum, for a UHQ 
microcavity pumped at 1450 nm. The pump-to-Raman 
conversion characteristics for first order Raman scatter- 
ing and the 2"'^ order Raman mode are shown in the 
inset. It can be seen that the first order mode does in- 
deed exhibit a square-root dependence on the launched 
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FIG. 10: Cascaded Raman Scattering in a 58-/im— diameter 
microcavity. The pump WGM is located at around 1450 nm. 
The insets show the pump-to-Raman conversion for first (left 
inset) and second order (right inset) Raman modes (measured 
on different microcavities) . The 1°* order Raman mode ex- 
hibits a square-root, and the 2"'' order Raman mode a linear 
pump-to-Raman conversion characteristic, in agreement with 
the theoretical prediction. Solid lines: A theoretical fit using 
equations 12. 
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FIG. 11: 5"* Order Cascaded Raman Scattering. The thresh- 
old of the cascades are plotted as a function of Stokes or- 
der (N). The data was acquired on a single microsphere of 
50-/xm— diameter excited using a WGM at 980-nm wave- 
length. The inset shows the efficiency of the cascades as a 
function of Cascading order, measured on several devices.. 



pump power. The solid line is a fit using equation 12. 
The higher order Raman mode exhibits the expected lin- 
ear increase with pump power. 

To study cascaded Raman scattering beyond 2"^^ order, 
experiments using a 980 nm wavelength pump [32] were 
employed. The shorter wavelength pump allowed the 
observation of up-to S'^'-order cascades (from 980-1300 
nm) owing to the reduced mode volume at shorter wave- 
lengths, and the higher Raman gain coefficient {gn cx 



1/A) . With less than 900 /i- Watts of launched fiber power 
up to fifth order Stimulated Raman Scattering was ob- 
served [32] , and the threshold and efficiency of the cas- 
cades measured. Figure 11 shows the Cascaded Raman 
threshold for a microsphere as a function of the order of 
the cascade, for fixed coupling condition. The solid line 
is a cubic fit (as predicted by equation 11), which yields 
good qualitative agreement with the experimentally mea- 
sured thresholds. The inset of figure 11 shows a double 
logarithmic plot of the measured differential conversion 
efficiency of the cascaded Raman scattering process. The 
efficiency of the cascades decreases as a function of Stokes 
order, as is theoretically predicted by equation 13. 



RAMAN SCATTERING IN ULTRA-HIGH-Q 
TOROID MICROCAVITIES ON-A-CHIP. 



Toroid microcavities allow on-chip integration of ultra- 
high-Q performance with other optical, mechanical and 
electrical functionality. In addition the wafer-scale fab- 
rication process allows precise dimensional control and 
parallelism. Figure 12 shows a scanning electron mi- 
croscope of a toroid microcavity. In addition to the 
design freedom and control brought about by micro- 
fabrication and integration possibilities, toroid microcav- 
ities also exhibit significant advantages in terms of their 
Raman-emission properties. Whereas microsphere cavi- 
ties show low threshold operation, their emission is in- 
herently multi-mode due to the presence of azimuthal 
modes, as shown in figure 7 of the previous section. In 
contrast, toroid microcavities on-a-chip have a signifi- 
cantly reduced mode spectrum, such that single mode 
Raman lasing can be observed. The reduced number of 
azimuthal modes is due to the cavity geometry, which 
supports only a few modes in the azimuthal direction 
due to the toroidal confinement [26]. Figure 13 shows the 
mode spectrum of a toroid microcavity. Only two higher 
order modes are present in the spectrum, and successive 
modes are separated by the free spectral range of the 
cavity, which in the figure is 10 nm. This is in contrast 
to a microsphere of identical principal diameter which 
would support 2i + I azimuthal modes weakly split by 
eccentricity. 

The suppression of azimuthal modes, has important 
consequences on the spectral emission properties of sin- 
gle and cascaded Raman scattering. Most notably, sin- 
gle transverse mode Raman oscillation can be observed 
over a large range of pump powers. Figure 14 shows 
a toroid microcavity Raman laser, which exhibits sin- 
gle mode oscillation at emission power levels up to 160 
fiWatts with high efficiency (45% at the critical point). 
The single-mode emission property is a significant advan- 
tage in practical applications of nonlinear optical oscil- 
lators. Toroidal and spherical microcavities having sim- 
ilar Q and outer diameter will differ in modal volume. 
Toroids, owing to the added transverse confinement will 
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FIG. 12: Scanning-electron-micrograph of a toroid microcav- 
ity. The principal diameter (D) and minor toroid diameter(d) 
are shown. 
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FIG. 14: Single mode Raman emission of a ca. 60- 
/im— diameter toroid microcavity.The pump is located at 1550 
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FIG. 13: The recorded mode spectrum (both transmission 
and reflection) of a 50-/im— diameter toroid microcavity, ex- 
hibiting a strongly reduced mode spectrum. The reduced den- 
sity of modes is due to the encreased azimuthal confinement, 
and is similar to spectra obtained with spheroids [26]. 



have a reduced mode volume and hence a lower Raman 
threshold. The degree to which the threshold is reduced 
will depend on the "aspect ratio" of the toroid or D/d 
where D is the outer or principal diameter and d is the 
minor diameter (see figure 12). The modal volume of the 
fundamental toroid mode in a D= 50 fxm toroid plot- 
ted versus d is provided in figure 15. The calculation 
uses a finite element mode solver as there are no analyt- 
ical expressions available for toroid mode volume. The 
case of d = 50 fim corresponds to a sphere. Overall, the 
toroidal mode volume can be seen to be lower than the 
sphere. There are also two distinct regimes of mode vol- 
ume behavior. In the first, the mode volume reduces very 
slowly as d is reduced {V cx d^^^). This regime features 
a weak lateral confinement very similar in nature to that 
of the sphere and related to the gentle transverse curva- 
ture of the toroidal dielectric boundary. In the second 
regime, the mode volume reduces very quickly as d is re- 
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FIG. 15: Numerically calculated effective mode volume of 
a toroid microcavity as a function of the cross-sectional di- 
ameter of the toroid (and fixed principal toroid diameter 
D=50/im). The circles denote the numerically calculated ef- 
fective mode volume for a fundamental toroid WGM in the 
1550-nm band for (TE) polarization. The inset shows the 
modeled TE-mode intensity distribution for a toroid with 2 
fim and 50 toroid-cross sectional diameter (d) using the 
finite-elements method. 



duced. This regime is characterized by strong, lateral, 
index confinement resulting from the toroidal boundary 
being comparable in diameter to the mode field. The 
inset to this figure provides mode fields for the toroid 
in these two regimes to further illustrate this idea. It is 
clear that substantial reductions in threshold are possible 
using these devices. We are currently characterizing the 
reductions possible experimentally. 
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VI. SUMMARY 

In summary we have experimentally and theoretically 
analyzed Raman oscillation in fibcr-tapcr-couplcd micro- 
spheres and microtoroids on-a-chip. A theoretical anal- 
ysis was presented using the coupled mode equations for 
the pump and Raman WGMs. Using these equations, 
the threshold condition for stimulated Raman scatter- 
ing was derived and the relative importance of waveg- 
uide coupling strength, mode volume and intrinsic res- 
onator Q were described. These theoretical dependences 
wore verified experimentally. Furthermore the analysis 
was extended to the case of cascaded Raman oscilla- 
tion and threshold and efficiency expressions were de- 
rived for higher-order Raman fields. This analysis re- 
vealed that odd and even order Raman lines exhibit dif- 
ferent pump-to-Raman emission characteristics. Even or- 
der Stokes fields are found to exhibit a linear increase in 
generated Raman power as a function of pump power, 
whereas odd-order Stokes fields exhibit a square root de- 
pendence. Analysis showed and experiment confirmed 
that the threshold for N-th-order cascaded Raman oscil- 
lation exhibits a cubic dependence on order and that the 
associated efficiency of the process scales inverse quadrat- 
ically with order. 

Microtoroids, in addition to having significant practical 
advantages with respect to their chip-based fabrication, 
have both spectral and power efficiency advantages in 



comparison to Raman oscillation in microspheres. Their 
stronger lateral confinement provides two distinct bene- 
fits. First, a drastic reduction in the complexity of the 
mode spectrum enabling single-mode oscillation in the 
microtoroid based device. Second, a controllable and re- 
duced mode volume so that for comparable Q factors, 
microtoroid devices should have lower threshold pump 
powers. We are currently investigating the latter advan- 
tage experimentally. 

The importance of fiber taper coupling in these mea- 
surements cannot be over emphasized. These inher- 
ently fiber-compatible waveguides provide exceptional 
coupling efficiencies to and from the ultra-high- Q de- 
vices. They are also indispensable in coupling to the 
microtoroid devices which reside near the surface of a sil- 
icon wafer. Using taper coupling, the lowest threshold 
observed in this study was 62 /x Watts of launched power, 
a value which is nearly 3-orders of magnitude lower than 
for free-space illumination of micro-droplets. 
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